A luminous fast radio burst that probes the Universe at 
redshift 1 


S. D. Ryder, K. W. Bannister?, S. Bhandari^?, A. T. Deller®, R. D. Ekers*’, 
M. Glowacki’, A. C. Gordon?, K. Gourdji?, C. W. James’, 
C. D. Kilpatrick??, W. Lu!?:!!, L. Marnoch! ^*^, V, A. Moss’, 
J. X. Prochaska!?: 4, H. Qiu’, E. M. Sadler!®?, S. Simhat’, 
M. W. Sammons’, D. R. Scott’, N. Tejos!”, R. M. Shannon?* 
!School of Mathematical and Physical Sciences, Macquarie University, NSW 2109, Australia 
? Astrophysics and Space Technologies Research Centre, Macquarie University, Sydney, NSW 2109, Australia 
3 Australia Telescope National Facility, Commonwealth Science and Industrial Research Organisation (CSIRO), 
Space and Astronomy, PO Box 76, Epping, NSW 1710, Australia 
“The Netherlands Institute for Radio Astronomy (ASTRON), 7991 PD Dwingeloo, The Netherlands 
5Joint institute for Very Long Baseline Interferometry in Europe, 7991 PD Dwingeloo, The Netherlands 
Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, VIC 3122, Australia 
"International Centre for Radio Astronomy Research, Curtin Institute of Radio Astronomy, 
Curtin University, Perth, Western Australia, Australia. 
8Center for Interdisciplinary Exploration and Research in Astrophysics, Northwestern University, 
Evanston, IL 60208, USA 
? Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208, USA 
10Department of Astronomy University of California, Berkeley, CA 94720, USA 
11 Theoretical Astrophysics Center, University of California, Berkeley, CA 94720, USA 
12 Australian Research Council Centre of Excellence for All-Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia 
13 Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA 
14 Kavli Institute for the Physics and Mathematics of the Universe, Kashiwa, 277-8583, Japan 
I5SKA Observatory, Jodrell Bank, Lower Withington, Macclesfield, SK11 9FT, UK 
16Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 
l"Instituto de Física, Pontificia Universidad Católica de Valparaíso, Casilla 4059, Valparaíso, Chile 


"Corresponding author: E-mail: rshannon@ swin.edu.au. 


Fast radio bursts (FRBs) are millisecond-duration pulses of radio emission 
originating from extragalactic distances. Radio dispersion on each burst is im- 
parted by intervening plasma mostly located in the intergalactic medium. We 
observe a burst, FRB 20220610A, in a morphologically complex host galaxy 
system at redshift 1.016 + 0.002. The burst redshift and dispersion are con- 
sistent with passage through a substantial column of material from the inter- 
galactic medium. The burst shows evidence for passage through additional 
turbulent magnetized plasma, potentially associated with the host galaxy. We 


use the burst energy of 2 x 10? erg, to revise the maximum energy of an FRB. 


Fast radio bursts (FRBs, (1, 2)) are transient radio sources that last a few milliseconds emit- 
ted by extragalactic sources. Free electrons along the path between the FRB source and the 
Earth impart a frequency-dependent time delay (dispersion) on the radio signal. This dispersion 
can be used to measure the column density of free electrons (quantified by the dispersion mea- 
sure, DM) between the FRB source and observer. FRBs localized to host galaxies at different 
redshifts exhibit a positive correlation between extragalactic DM and host redshift, known as 
the Macquart relation (3). This relation has been used to measure the cosmic baryon fraction 
and the expansion rate of the Universe (4). This relation has been measured using identified 
FRB host galaxies at relatively low redshifts, z < 0.5. Some unlocalized FRBs (with unknown 
host galaxies) have DMs consistent with z > 1 (5); however an FRB associated with a galaxy 
at z = 0.241 had a high DM that would have implied z > 1 (6). This indicates that estimates 
of redshift from DM alone can be misled by plasma within the host galaxy, which also imparts 


a contribution to the DM, in addition to that of the intergalactic medium. 


Observations of FRB 20220610A 

FRBs have been searched for (7) and localized (8) using the Australian Square Kilometre 
Array Pathfinder (ASKAP, (9)), a radio interferometer in Western Australia comprising 36, 12- 
m antennas. Each antenna is equipped with phased-array receiving systems, which provide 36 
beams across the focal plane, covering approximately 30 square degrees. FRBs are searched for 
in real time using the incoherent sum of the intensities of each antenna in each beam. When an 
FRB is detected, voltage buffers are downloaded, correlated, calibrated, and imaged, enabling 
the position of the burst to be measured to an absolute precision of typically a few tenths of an 
arcsecond (10, 11). 

We detected FRB 20220610A in observations around the previously known burst FRB 20220501C, 
but the two bursts are not related (77). The observations were centred at a frequency of 
1271.5 MHz and had a time resolution of 1.18 ms. The dispersion of FRB 202206104 indi- 
cated a DM of 1458.15*02? pc cem^?. This is higher than all but one of the 55 FRBs previously 
observed using ASKAP (4). A dedispersed dynamic spectrum of the burst is shown in Figure 1, 
and properties of the burst are listed in Table 1. The burst does not show the 10-100 MHz modu- 
lation in the spectrum characteristic of many lower DM, high Galactic latitude ASKAP-detected 
FRBs (7). 

The 25 dispersive sweep of the burst across the instrument bandwidth and 2.4 s latency in 
the detection system resulted in only the lowest 88 MHz of the burst being captured in the 3.1 s- 
duration voltage buffer. This was sufficient to localize the burst to a precision of 0.5 arcsec. We 
used the voltage data to reconstruct the high time resolution and polarimetric properties of the 
burst (//). After correcting for dispersive smearing, the burst shows an exponentially decreasing 
tail (Fig. 1D), which is consistent with scatter broadening due to turbulence in intervening 
plasma (72). We measure the pulse broadening time to be 0.511 + 0.012 ms at a reference 


frequency of 1147.5 MHz, assuming a v 4 frequency dependence. 


Ordered magnetic fields in astrophysical plasmas add additional, polarization-dependent 
dispersion. This manifests as wavelength-dependent variation in the linear polarization position 
angle, referred to as Faraday rotation (13). The burst exhibits Faraday rotation, with a rotation 
measure (RM) of 215 + 2rad m ?. After correcting for this Faraday rotation, we find the burst 
had a linear polarization fraction of 963-1 %. The high fractional linear polarization allows us to 
place a 6796 upper limit on the Faraday dispersion (og < 0.6 rad m^?) induced by fluctuations 
in rotation measure in intervening turbulent plasma. Higher levels of Faraday dispersion have 
been detected for other FRBs (74). The burst also shows modest fractional circular polarization 
of 1023-196. While instrumental artifacts can induce spurious circular polarization, we do not see 
any correlation between Stokes polarization parameters U and V in the spectrum, which would 
be expected for an instrumental effect (77). The FRB was located approximately 4 arcmin from 
the beam center, which makes off-axis leakage effects less likely (75). Circular polarization 
has been observed in some FRBs and could either be intrinsic to the burst (/6) or result from 


propagation through relativistic plasma in the immediate source environment (17). 


Host-galaxy properties 

We performed follow-up ground-based optical and infrared observations with the Very 
Large Telescope (VLT) and the W. M. Keck Observatory to identify and characterize the host 
galaxy of FRB 20220610A (7/1). The images (Fig. 2A-C) show an object coincident with the 
source that has an extended, multi-component morphology. We label the optical source that 
overlaps the radio position of the FRB as component (a), and two adjacent sources as compo- 
nents (b) and (c) (Fig. 2A). We use a Bayesian method to assess the chance of coincidence 
between transients and host galaxies (/8), finding greater than 99.99% confidence that the FRB 
is associated with component (a). 


We performed broad band optical and infrared spectroscopy of components (a), (b) and (c) 


(Fig. 2D-E) (11). We identify two emission lines in the spectra as the [O 1] 3726 and 3729 A 
doublet, most prominently in component (b) (Fig. S2). From this we measure the redshift of 
each component, finding they are all consistent with z = 1.016 + 0.002. 

We estimate the total mass of the three components combined to be 10!? solar masses, with 
a star formation rate of ~ 0.42 solar masses per year (/7). These values, in addition to the host 
metallicity and star formation history are consistent with those of nearby FRB hosts (/9, 20), but 
the source morphology is markedly different. Properties of the host galaxy are listed in Table 1. 

The presence of two bright components (a) and (c) separated by 2.0 arcseconds (which 
corresponds to a distance of 16 kpc at that redshift), and the diffuse feature (b) between them, 
is consistent with two galaxies interacting or merging, or a compact galaxy group. It is also 
possible that the morphology is due to internal structure within a single galaxy; at these redshifts 
about half of all galaxies have clumpy morphologies (2/). We regard the latter possibility as 
unlikely, due to the large spatial separation between the components. Only component (a) is 
detected in the near-infrared (J€,-band) image (Fig. 2C), indicating it hosts an older stellar 
population than the other components. Component (a) is also displaced from the centroid of the 
total optical light in g- and R-bands, contrary to what would be expected if it was the nuclear 


bulge of a single galaxy. 


Extending the Macquart Relation 

We used the measured properties of FRB 20220610A to probe the Macquart relation to 
z ~ 1, by comparing with predictions for its DM based on previous fits to the relation at 
z < 0.522. Figure 3 shows the relationship between DM and redshift for the FRBs detected by 
ASKAP (4). We restrict our analysis to the ASKAP sample (4) to minimize observing-system- 
dependent selection effects. We do not re-fit the Macquart relation. Doing so in an unbiased 


way would require analyzing the entire updated FRB sample from ASKAP. 


After subtracting a model Milky Way foreground contribution to the observed DM (11), 
we estimate the non-Galactic DM of FRB 20220610A to be ~ 1376 pc cm ^, indicating a high 
column density of ionized gas between the FRB and Earth. This is higher than the DM ex- 


pected from the Macquart relation, by approximately 650 pc cm? 


, and is a 2.40 excess (11). 
If the excess is real and originates from the host galaxy, the implied electron column density is 
13007355 pc cm^? in the host rest frame (77), with the uncertainty reflecting the intrinsic varia- 


tion in the contribution from the intergalactic medium. 


Interpretation of the Dispersion Measure 

We use the scatter broadening and Faraday rotation of FRB 20220610A to investigate the 
properties of plasma at z ~ 1. The possible dispersion excess could arise from any combination 
of gas in the immediate vicinity of the source, the interstellar medium of the host galaxy, or fore- 
ground gas along the line of sight, each of which could potentially harbor turbulent magnetized 
plasma. The absence of strong (> 104 rad m^?) Faraday rotation or detectable depolarization is 
unlike the repeating FRB 20190520B, which also shows excess dispersion (6, /4). We suggest 
the excess dispersion for FRB 20220610A originates in a less magnetoionically active plasma 
than FRB 20190520B, such as in the interstellar medium of the host galaxy, rather than in the 
circumburst media hypothesized for other sources (77). Models of galaxy interstellar media im- 
ply that the DM of a typical spiral galaxy is unlikely to exceed a few hundred pc cm ? except in 
edge-on systems (22, 23). Higher DM values can plausibly be produced by high-density clumps 
of gas within the host galaxy, particularly at z ~ 1 where galaxies have a substantially higher 
fraction of their baryons in gas (rather than stars) compared to z ~ 0 (24). Alternatively the 
dispersion could originate from structure in the foreground intergalactic medium, or additional 
ionized material associated with the possible galaxy merger between components (a), (b), and 


(c). 


Our DM analysis confirms inferences from other techniques (25) that the gas of the inter- 
galactic medium is highly ionized. The detection of an FRB at z > 1 allows us to study the 
ionized plasma towards, around, and within the host galaxy. We expect a sightline to z = 1 to 
intersect the halos of several galaxies similar in mass and size to the Milky Way (26), yielding 


further insight into their properties (27). 


Interpretation of the burst energy 

The measured bandwidth-averaged fluence of FRB 20220610A is 45 + 5Jyms (1 Jy = 
10 73 eres ! cm ? Hz~'), implying an isotropic-equivalent spectral energy density of (6.4 + 0.7) x 
10°” erg Hz~! and a burst energy of (6.4 + 0.7) x 10‘ erg (11). We derived the burst energy 
by assuming an intrinsic 1 GHz bandwidth, and did not apply a redshift-dependent correction 
to the burst spectral energy distribution (i.e., K-correction) (77). This value exceeds the char- 
acteristic maximum energy Emax derived by previous FRB population models by a factor of 
3.5 (11, 28). It is unknown whether FRBs are emitted isotropically or only over a limited beam- 
ing angle, which would affect the inferred energetics and could vary between FRB sources or 
repeat FRBs from single sources. Assuming isotropic emission, we have re-fitted the FRB burst 


energy distribution by adding FRB 20220610A to the sample of FRBs used by (28), finding the 


best-fit Emax to increase by a factor of 2.7, to 10417 + 9? erg (i.e., the maximum energy density 
becomes logy, E, [erg Hz !] = 32.7 + 0.2). In Figure 4, we contrast the burst to the fluence 
of the brightest radio pulse observed from a galactic magnetar, which is a factor of 10? less 
luminous than the burst observed from FRB 202206104, and the wider sample of FRBs. 

If a K-correction is applied, assuming a spectral index for FRB emission similar to that of 
the burst population found by ASKAP (1/1, 29), we find the burst energy integrated over the 
instrument bandwidth is ~ 2 x 10? erg, which is higher than most localized FRBs (Figure 


4) (3, 4, 6, 8, 17, 27, 30—37). This constrains emission models of FRBs, because the electric 


field strength at the source can be estimated independently of the beaming angle. The cal- 
culation assumes no amplification by gravitational or plasma lensing. In one class of mod- 
els, FRB emission is produced near the surface of a neutron star. From the luminosity of the 
burst ~ 3 x 10*° ergs! in the host galaxy's rest frame, we infer an electric field strength of 
4.2 x 10? (r/10km) ! V m^, for a linearly polarized wave, where r (~10 km) is the curvature 
radius of the neutron star's magnetic field. This is a few percent of the Schwinger critical field 
strength at a neutron star surface, in which an electric field aligned parallel to the local magnetic 
field, would be quickly screened by copious electron-position pair production (38). This would 
suppress the FRB rate above the Schwinger luminosity of approximately 2 x 10?" erg s^! (38). 
In another class of models, FRBs are produced in a shock driven by relativistic ejecta associ- 
ated with the flare of a highly magnetized neutron star interacting with the neutron star wind. 
In these models, the radiative efficiency in the shock is very low (< 10^?) and hence the re- 
quired energy in the ejecta would be > 10?" erg, with the total flare energy being even higher. 


FRB 202206104 and other high luminosity FRBs challenge both model classes. 
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Table 1. Properties of FRB 20220610A and its host galaxy. The fluence was derived from 
filterbank data that includes the full ASKAP bandwidth. UTC is the Coordinated Universal 
Time standard. S/N is the signal to noise ratio. FWHM is full width at half maximum. The 
brightness of the host is expressed in magnitudes (mag) on the AB system. The star formation 
rate (SFR) is in solar masses per year (Mo yr~'). The integrated star formation in the last 100 
Myr is listed in the row 100 Myr SFR. The metallicity of the host galaxy, Z, expressed as a 
(logarithmic) fraction of Solar metallicity Zo is listed in the row log(Metallicity). 


Measured burst properties 


Dispersion Measure (DM) 1458.15 77 pean? 
Topocentric arrival time at 1104 MHz (UTC) — 2022-06-10 22:26:44.313 
Fluence 45 + 5 Jy ms 
Right Ascension (J2000 equinox) 23^24"17.5695 + 0.0405 
Declination (J2000 equinox) —33°30'49".37 + 0".50 
Galactic Longitude 8.83954? 
Galactic Latitude — (0.18569? 
Incoherent Detection S/N (1104 — 1440 MHz) 29.8 
Image S/N (1104-1152 MHz) 81 
Rotation Measure 215 + 2 rad m~? 
Inferred burst properties 
Intrinsic width (FWHM) 0.41 + 0.01 ms 
Implied FRB isotropic energy density 6 x 10? erg Hz! 
Milky-Way disk DM contribution 31 pc cm~’ 
Measured host galaxy properties 
Redshift 1.016 + 0.002 
Photometry 
Band (Central wavelength, Å) Magnitude (mag) 
g (4700) 24.15 4 0.07 
V (5510) 23.89 + 0.13 
R (6580) 23.78 + 0.06 
I (8060) 22.17 + 0.07 
z (9620) 21.95 + 0.12 
J (12200) 21.97 + 0.07 
Ks (21,460) 22.08 + 0.12 
Inferred host galaxy properties 
Mass-weighted age 1.024185 Gyr 
log(Stellar mass) 9.987557 Mo 
log(Total mass) 10.11848 Mo 
100 Myr SFR 0.42+9:83 Mo yr 
log(Metallicity) —0.11*947 Z/Zo 
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Figure 1. Radio observations of FRB 20220610A. (A and B): The burst as observed in the 
incoherently summed data stream used by the real-time detection system. The pulse width is 
dominated by intra-channel dispersion smearing. Panel A shows the integrated flux density (S) 
of the pulse profile as a function of time (t) while Panel B shows the dedispersed burst dynamic 
spectrum as a function of frequency (Freq). The vertical dashed line in panel B indicates the 
frequency range (<1200 MHz) of high time resolution data that were saved by the pipeline; 
data at higher frequencies were lost due to latency issues (see text). (C and D): High time 
resolution data (frequencies <1200 MHz) produced from the raw electromagnetic field samples 
saved from each telescope after the real-time trigger (77). Panel C shows linear polarization 
position angle V during the pulse. Panel D shows the integrated pulse profile in total intensity 
I (black), linear polarization Leg (red), Stokes V (blue), and a maximum likelihood model of 
Stokes I, M; (orange). The residuals Rz (pink) between J and M; are shown offset from zero 
by —20 Jy for clarity. The polarisation position angle V is defined such that linear polarization 


is in one Stokes parameter and hence V z 0. 
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Figure 2. Optical and infrared observations of the host galaxy of FRB 20220610A. (A) 


7505 


g-band VLT image. White ellipses outline the apertures used for the photometry of each com- 
ponent. The larger unlabeled white circle is the aperture used for the entire combined system. 


The FRB localization and uncertainty (6896 confidence) are indicated by the black ellipse in 
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panels A-C. (B) R-band image, with the slit locations used for VLT spectroscopy superimposed 
in white. (C) Ks band image. The same scale bar applies to panels A-C and shows the angular 
scale and corresponding projected physical scale at the measured host redshift. The color bar 
below each panel indicates relative count rates in each processed image. (D) Two-dimensional 
VLT spectroscopy at a slit position angle on the sky of 45° covering components (a) and (b) 
as marked in Panels A and B. White arrows indicate the two lines due to [O II] with rest wave- 
lengths of 3726 and 3729 A. (E) One-dimensional spectrum (black solid line) and its uncertainty 
(red dashed line) extracted from the image shown in Panel D, centred on the peak of the [O 11] 


lines using an aperture width of 1.8". The blue line indicates zero relative flux in arbitrary units. 


21 


3000 


E 
B 
oO 
a 
1800 E 
o 

o 
ui c 
= S 
a o 
1200 2 
2 
S 
o 
o 
600 = 


Figure 3. Relationship between redshift and extragalactic dispersion measure for FRBs. 
Data are from the ASKAP incoherent sum survey (4). The extragalactic dispersion (DMgo) is 
the contribution to DM after subtraction of a model for the Milky Way (//). The localized FRBs 
are shown as red circles at the host galaxy redshifts. The color scale indicates an estimated 
detection probability assuming an increased maximum energy density log;; E, [erg Hz !| = 
32.7 + 0.2 (found in this work). White contours show 50% (dotted), 90% (dash-dot) , and 
99% (dashed) of the probability. FRB 20220610A is plotted with the white star. The repeating 
FRB 20190520B, shown as the black cross, was not used in the population inference due to 
varying selection effects between the discovery instruments (28). The uncertainties for the 


measurements are smaller than the symbol sizes. 
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Figure 4. Logarithmic plot of fluence as a function of redshift for localized FRBs. The 
magenta star is FRB 20220610A. Other notable FRBs are labeled using their Transient Name 
Server designations, omitting the FRB prefix for brevity. The curved solid and dashed con- 
tours indicate the energy density in units of erg Hz !. The dash-dotted lines show the detection 
sensitivity of the current ASKAP incoherent sum FRB search system, and Five-hundred-meter 
Aperture Spherical Telescope and Square Kilometre Array telescopes (labeled FAST/SKA). Re- 
peating FRBs are shown at multiple points at the same redshift but different fluences. Plotting 
symbols are indicated in the legend. FRB localization instruments are noted in the top right; 
DSA is the Deep Synoptic Array. The green star shows the expected fluence of the FRB-like 
burst from Galactic magnetar soft gamma repeater (SGR) 1935+21 (39, 40), if it were placed at 


the distance of the host galaxy of FRB 20180916A. The repeating FRB 20201124A is plotted 
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twice for observations from ASKAP (labeled “20201124A (A)") (77) and the 25-m Stockert 
radio telescope (labeled “20201124A (S)’’) (33). We also show the initial detection of the re- 
peating FRB 201211024 (labeled “20121102A (D)") (41), and a sample of follow up bursts de- 
tected at Arecibo (labeled “20121102A (R)’”) (35); and bursts from two low-redshift repeating 
FRB sources: FRB 20180916A (36) and FRB 20200120E (37). Burst properties are reported 
from ASKAP (3,4, 8,27); DSA-10 (30) and FAST detections of repeating FRB 20180301A (31) 
and FRB 20190520D (6), in addition to the initial detection of FRB 20180301A with the Parkes 
radio telescope (labeled “20180301A (D)’) (32). There are biases in FRB searches which may 
result in the first detections of repeating sources having different properties than subsequent 


detections, particularly those made with other facilities (42). 
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S1 Materials and Methods 
S1.1 Radio observations 


FRB 20220610A was detected as part of the Commensal Real-time ASKAP Fast Transient 
(CRAFT) survey (43), which uses ASKAP (9) to search for dispersed radio transients such 
as FRBs. The PAF receivers on ASKAP enable 36 beams to be digitally synthesized across the 
focal plane of the antenna. The FRB detection system on ASKAP at the time of observation 
searches intensity data, incoherently summed across telescopes, in real time. Due to constraints 
on performance of the search system, only a subset of the 36 antennas were included in the real- 
time sum, with trade offs between the number of antennas included, and the time resolution of 
the data searched. During the observation in which FRB 20220610A was detected, the searches 
were conducted using 22 of the antennas, with a time resolution of 1.182 ms. In all cases a 
total bandwidth of 336 MHz was searched with a spectral resolution of 1 MHz. The search data 
streams were recorded to disk in pulsar filterbank (44) format for subsequent analysis whenever 
a burst was detected. Details of the search system have been described previously (8). 

FRB searches with ASKAP are intended to be commensal with other observations. How- 
ever some observations are also conducted in a mode termed filler, during times where standard 
observations cannot be scheduled. The observations reported here were conducted at a central 
frequency of 1271.5 MHz with the configurable beam in a hexagonal close packed configura- 
tion and beam centers separated by 0.9 degrees. FRB 20220610A was detected during filler 
observations while pointed towards FRB 20220501C, a previously-detected FRB (45). The 
different sky locations (the burst positions, with precisions better than one arsecond, are sep- 


3 and 


arated by approximately 1.5 deg) and burst dispersion measures — 449.5 + 0.2 pc cm7 
1458.15*02? pc cm^?, respectively — rule out the two bursts being from the same repeating 


Source. 


The burst was identified with a S/N of 29.8 in the primary beam. The ASKAP beams 
oversample the focal plane, so the burst also appeared at S/N> 5 in four adjacent beams. The 
burst width, at full width half maximum as measured in the filterbank, was 5.6 + 0.4ms. This 


is consistent with the 6.0 ms burst width expected from intrachannel dispersion smearing (46). 


S1.2 FRB localization 


From each ASKAP antenna, 3.1 s of voltage data was saved after being triggered by the real- 
time detection system. Due to the dispersive sweep of 2.04 s across the 336 MHz-wide ASKAP 
band, combined with the 2.4s latency of the detection system, much of the FRB had already 
been lost from the voltage buffer before it was frozen and downloaded. Below a frequency 
of 1192 MHz, the FRB signal remained present in the saved data. Approximately 2.5 h after 
the FRB detection, a voltage download was triggered while observing the bright radio source 
PKS 0407-658, to be used for calibration purposes, followed shortly after by a voltage download 
on the Vela pulsar (used as a polarization calibrator). 

Our localisation of FRB 20220610A used the ASKAP FRB astrometric pipeline (8, 47). We 
correlated a small time window matching the dispersed FRB itself (referred to as the gated 
dataset), as well as the full 3.1s time window from the FRB voltage download (the field dataset) 
and the voltages downloaded while observing PKS 0407-658 (the calibrator dataset). The volt- 
ages downloaded while observing the Vela pulsar were also correlated (the polcal dataset). The 
downloaded voltage data were processed using the CELEBI pipeline (48), which is based on the 
same procedures as for earlier ASKAP FRBs (6, 47). Voltage data were correlated using the 
DIFX software correlator (49), to produce the four visibility datasets. 

Calibration was derived from the calibrator dataset and applied to the gated, field, and polcal 
datasets following previously described procedures (47). None of the 22 antennas had to be 


removed from analysis for reasons such as irregular visibility amplitudes that occasionally affect 


a subset of downloaded CRAFT data. After calibration, each dataset was imaged using CASA 
(50). Portions of the spectrum affected by radio frequency interference (for the field and polcal 
datasets) and the lack of FRB emission due to the limited voltage buffer size and FRB dispersion 
(for the target dataset) were flagged and discarded prior to imaging. The intrinsic width and 
scattering time of the FRB were fitted using a Bayesian method described below (section|S1.4). 
Positions and corresponding uncertainties were measured in the image plane using the AIPS 


task JMFIT (51). 


The target position was measured from the FRB image to be right ascension 23524" 175.5594 
05.006, declination —33^30'49".87 +0”.07, J2000 equinox. We used an ensemble of seven back- 
ground radio continuum sources from our 3.1 s field image to estimate any systematic shift in 
the FRB position, following published methods (52). These seven sources were selected be- 
cause they are all the detections in this field image that were unresolved in the National Radio 
Astronomy Observatory Very Large Array Sky Survey (NVSS) (53). All these sources were 
also present in the Rapid ASKAP Continuum Survey (RACS) catalog (54). Neither catalogue is 
ideal for this comparison; RACS is at a lower frequency than our observations and the catalogue 
used in this analysis is affected by astrometric errors of up to <1 arcsecond (54) while NVSS 
is at higher frequency and lower angular resolution. We measure a small astrometric correction 
from each catalogue, which are consistent at the 1c level: —0.16" 0.5" and 0.85" +0.5” in right 
ascension and declination respectively for the RACS catalogue positions, and 0.39" + 0.5" and 
0.14" 4 0.5" in right ascension and declination respectively from the NVSS catalogue positions. 
This offset is small and consistent with previous results (52). 

Because we have no reason to prefer one catalogue over the other, we take an unweighted 
average of the two offsets as our final correction: 0.12" and 0.50" in right ascension and dec- 
lination respectively, with an uncertainty on the correction of 0.5" in each coordinate. This 


is equivalent to using an unweighted average of the RACS and NVSS reference positions for 


each source when performing the comparison to the field image positions, without reducing 
the corresponding reference position uncertainty. The uncertainty on our systematic position 
correction dominates the final position uncertainty. After applying this astrometric correction, 


the final FRB position is 23524™17°.569 + 05.040, —33?30'49".37 + 0".50 (J2000). 


$1.3 Beamforming and spectral-temporal polarimetry 


The voltage data were also used to reconstruct a complex-valued time series of the electric 
field in two linear polarizations, from which we derived the Stokes 7, Q, U, V parameters and 
measurements of the polarization properties of FRB 202206104 at a time resolution of 3 ns. The 
beamforming process used to achieve this has been previously outlined (76, 48). This pipeline 
(CELEBI, (48)) applies per-antenna time delays to account for the difference in signal arrival 
times and uses a polyphase filterbank inversion to obtain a single complex, high frequency 
resolution spectrum for each polarization per antenna. Calibration solutions were then applied, 
before spectra were summed across all antennas and derippled. Differential gain and phase 
between the two linearly polarized feeds was derived using the polcal dataset of the Vela pulsar 
(PSR J0835—4510). The FRB data were then coherently dedispersed (using an initial DM of 
1457.6 pc cm?) and a Fourier transform inversion applied. Finally, the Stokes products were 
calculated and the data averaged in time and frequency to 1 us resolution. 

To find the S/N-maximising DM, the data were coherently dedispersed to several trial DMs 
about the nominal value of 1458 pc cm ?. For each trial DM, Stokes J was averaged at 10 us 
resolution, and the S/N calculated for sliding boxcar windows of varying width. The peak value 
of S/N was found at 1458.15* 02? pc cm^?, with the uncertainties — consistent with our initial 


estimate of 1457.6 pc cm ? 


— corresponding to the range producing a S/N within 1 of the 
peak. These uncertainties reflect the intrinsic burst structure and not random noise: doubling 


the amount of noise in the data (by adding random Gaussian deviates with standard deviation 


equal to that found in off-pulse data), and re-performing the procedure 1000 times, results in 


68% of the data falling in the range of 1458.154995 pc cm. 


S1.4 Burst scatter broadening 


We use a nested sampling approach (55), implemented using the BILBY (56) Python library to 
determine the maximum a-posteriori model of the pulse profile. Assuming uniform priors on 
the pulse broadening time 7 and full width at half maximum w, and using a uniform weighting 
across the burst profile while assuming Gaussian errors to calculate and maximize the posterior 
probability of the model, the peak posterior distribution parameters are 7 = 0.511 + 0.012 ms 
and w = 0.407 + 0.012 ms. The NE2001 model of the Galactic electron density variations 
(57) predicts the Galactic scatter broadening, Tmw, to be 48.3 ns at 1.1475 GHz, so we infer 
the scatter broadening is dominated by the host environment, consistent with a large host DM 
contribution. Fitting our model of pulse broadening to four 22 MHz sub-bands and assuming a 
T x v ^ scaling yields a consistent result for the scattering time. 

Taking the measured redshift of 1.016 along with our reference frequency of 1.1475 GHz, 
and assuming a scattering index of 7 x v ^, we scale to the host galaxy rest frame (indicated 
with ^) to find Ti 454, = 1.88 + 0.02 ms, and w’ = 0.202 + 0.006 ms. In Table[S1]we compare 
our measurements to FRB 20190520B, another FRB with large excess DM above the Macquart 


relation, with the same redshift corrections applied. 


S1.5 Burst rotation measure 


After polarimetric calibration, we found that FRB 20220610A was highly linearly polarized, 
and that the Stokes Q and U parameters showed frequency-dependent variation that we attribute 
to Faraday rotation in cold plasma. Because the polarization of the burst is constant in time, we 


measure the rotation measure (RM) of the burst after averaging the de-dispersed burst in time 


to form one-dimensional (time-independent) spectra of the burst Stokes Q and U parameters: 
Q(v) and U (v) (58). We calculated these by integrating over a 1.9 ms gate around the peak of 
the pulse and used a Bayesian inference method described previously (6, 59). In brief it fits a 
model of the expected frequency-dependent variations in Stokes Q and U, while marginalizing 
over intensity variations in individual channels. We find the best-fitting rotation measure to be 
215 + 2 rad m ?. The pulse-averaged spectrum of the burst and the best-fitting model for the 
rotation measure-induced change in the polarization position angle V are presented in Figure 

The burst shows strong linear polarization with no evidence for polarization position angle 
changing across the pulse. The noise statistics of this quantity are Gaussian as it is a linear 
combination of components that also have Gaussian noise distributions, unlike the more gen- 
eral estimate of linear polarization L = JQ? +U? (47, 60), which requires a debiasing of 
linear polarization intensity (60). Linear and circular polarisation fractions were calculated by 
integrating (using a box car window) the polarized flux over the on-pulse region, and normal- 
izing by the total intensity. Uncertainties in linear and circular polarisation were determined by 
measuring the root mean square of the off-pulse baseline and scaling by the square root of the 


number of samples averaged in the spectrum. 


S1.6 Optical and infrared imaging 


We obtained images of the location of FRB 20220610A using the FOcal Reducer and low 
dispersion Spectrograph 2 (FORS2) on Unit Telescope 1 (UT1) of the European Southern 
Observatory (ESO) Very Large Telescope (VLT) on Cerro Paranal. On-chip binning of 2x2 
yielded a pixel scale of 0.25" over the 6.8’ x6.8’ field of view. Five exposures of 400 sec in the 
R_SPECIAL filter (61), with offsets of +10” between each, were obtained on 2022 July 1 UTC, 


while a total of 10 exposures of 600 sec in the g filter were obtained across 2022 Oct 2-4 UTC, 


in seeing of 0.9 to 1.0”. The FORS2 images were processed as described elsewhere (3), with the 
photometry calibrated against Dark Energy Camera Local Volume Exploration survey data (62) 
for the g-band images, and using the ESO FORS2 Quality Control level-1 database for R-band. 

Additional optical imaging of the host of FRB 20220610A was obtained with the Deep 
Imaging Multi-Object Spectrograph (DEIMOS) on the Keck-II telescope on Maunakea, Hawaii 
on 2022 September 24 UTC in the V, J, and z filters for 3x350 s in each. Atmospheric con- 
ditions were near-photometric with an average seeing of 0.9" and the field was observed at an 
average airmass of 1.7. All imaging was reduced following standard procedures in the Python- 
based package POTPYRI (63). Photometric zero points for the stacked frames were derived 
using Landolt photometric standards (64) transformed into the DEIMOS photometric bands. 

Near-infrared imaging was carried out with the High Acuity Wide-field K-band Imager 
(HAWK-I) camera (65) in combination with the ground-layer adaptive optics system on Unit 
Telescope 4 of the VLT, covering a 7.5'x7.5' field at 0.106"/pixel. Observations in the K, 
filter were taken on 2022 July 24 UTC, and in the J filter on 2022 September 29 UTC in 
clear conditions and delivered image quality of 0.4 to 0.5". For each filter 15 integrations of 
10 sec were averaged at each of 16 random offset positions within a 20" box. Image mosaics 
were processed using ESOREFLEX (66), with photometric calibration undertaken using the Two 
Micron All-Sky Survey point source catalogue (67). 

For the DEIMOS imaging, aperture photometry (reported in Table 1) was performed using a 
3" host radius and a 4" to 7.5" background annulus. For the VLT imaging two approaches were 
used to obtain the photometry reported in Table 1: first, a single circular aperture with a radius 
of 3" encompassing the entire source; and second, three ellipses defined by eye (Fig. 2A) from 
the R-band imaging, to individually measure components (a), (b) and (c). The results are listed 
in Table [S2] Both approaches, including background subtraction, used the SEP package (66). 
Table|S2]also lists the positions of the clumps as defined by the apertures. 


S1.7 Optical and infrared spectroscopy 


Spectroscopy covering the entire ultraviolet through near-infrared wavelength range of 0.3— 
2.48 uum was conducted with the X-shooter spectrograph (69) on Unit Telescope 3 of the VLT. 
Slit widths of 1.0", 0.9", and 0.9" yielded spectral resolving powers of 5400, 8900, and 5600 in 
the ultraviolet/blue (UVB), visible (VIS), and near infrared (NIR) arms, respectively. On-chip 
binning of 2x2 was employed in the UVB and VIS detectors due to the faintness of the target, 
without affecting the spectral resolution in these arms. 

Observations of the target were made at 2 slit position angles: the first on 2022 July 29 UTC 
at a position angle of 45° covering components (a) and (b); and the second on 2022 Septem- 
ber 23/24 UTC at a position angle of 90? covering components (a) and (c), as illustrated in 
Figure 2 (B). At each of these position angles, a total of 7200 and 7548 sec respectively of inte- 
gration time was accumulated in clear conditions and seeing 0.5 to 1.2", with frequent nodding 
of the target by +2.5” along the slit to assist with sky subtraction. Observations of the B9 V 
stars HD 138940 and HD 168852 were obtained to enable removal of telluric features, and of 
the spectrophotometric standard stars EG 21 and LTT 987 to allow relative flux calibration. The 
spectra were reduced with the PYPEIT data reduction pipeline (70) using standard techniques 
for image processing and wavelength calibration. The VIS and NIR science exposures were 
individually processed and the sky-subtracted 2D spectral images were co-added to maximize 
S/N. No emission lines were detected in the UVB arm spectra. Figure 2 (D-E) shows the 2D 
image and 1D extracted spectrum at the slit position angle of 45°. 

From the spectra of components (a), (b) and (c) we identified emission lines at 7510 and 
7515 A as the [O 11] 3726 and 3729 A forbidden transitions of ionized oxygen (Figure 2). Com- 
bining the measurements of these lines we determine the redshift of the host galaxy components 


(a)-(c) to be 1.016 + 0.002. 


S1.8 Dispersion-measure analysis 


Our DM analysis is based on a four-component model (5), 


DMhost 
14z' 


DM = DMww + DMhalo + DMeosmic F (S1) 


where DMyrw is the Milky Way disk DM contribution, DMyj4j, is the Milky Way halo contri- 
bution, DMeosmic 1s the cosmic DM contribution and DM, is the host DM contribution. 

We adopt a contribution from the Milky Way interstellar medium (ISM) of DMmw= 31 
pe em ^? (57) (cf. 13.6 pc cm ? predicted by the YMW model (71)), and a halo contribution 
of DMnato = 50 pc em ?, which is on the low side of predictions (26), but consistent with the 
observation of nearby, low-DM FRBs (72). The contribution from the intergalactic medium 
(IGM), including filaments and intervening halos, is a probabilistic distribution about the mean 
density of ionized gas in the universe (73). We adopt a previously described model (3), which 
is implemented in the FRB Python library (74). We use cosmology from the results of the 
Planck spacecraft (75), in particular the Hubble constant Ho = 67.4 kms ! Mpc !, physical 
baryon density parameter Nh? = 0.0224, and assume a feedback parameter F = 0.32. The 
parameter F governs the clumpiness of cosmic gas distributions, and hence the variance of 
DM.osmic about its expectation value. PF is the least well constrained parameter in the model, 
with log, F > —0.89 at 99.7% confidence (45). At the redshift of 1.016 + 0.002, this model 
predicts a mean intrinsic DM cosmic of 805+272 pc em ? (1o uncertainty). We attribute the re- 
maining contribution to DMhost, which includes the FRB host galaxy halo and ISM, and any 
plasma in the immediate vicinity of the FRB source. This contribution is suppressed by a factor 
of 14- z (here, ~ 2), because the dispersion is induced at higher frequencies than those observed. 

Using Eq. (SI), we estimate the probability distribution of DMnost for FRB 20220610A, 
with the results shown in Figure Because DM.og,;, is upwardly skewed, with the most 


likely DM being below the mean, and a small probability of very large excess DMs due to 
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intersections with cosmic structures, our estimate for DMhpost is left-skewed. The most likely 


3. with a 6896 confidence interval of 


value in the host rest-frame is DMhost = 1304 pc cm7 
(981, 1475) pc cm ?. However, we cannot exclude the possibility of a negligible host contribu- 
tion: 7% of the probability density lies in the unphysical region DMpost < 0. 

We also calculate the equivalent distribution for FRB 20190520B. This FRB originates from 
a z = 0.241 host galaxy (6). Its observed DM»919o52op = 1204.7 pc cm ^? implies a host galaxy 
contribution in the observer frame of 9U8* S. pc em ?. Using our model for DMeosmics DMhalo. 
and adopting DMyw = 60 pc cm ? (6), we find a rest-frame DM of 1172 pc cm ? with a 
68% confidence limit of (1104, 1197) pc cm ?. The predicted range is smaller due to the lower 
redshift of FRB 20190520B, and hence lower absolute variance in DMeosmic. This agrees with 
previous interpretations (76) that the DM of this FRB is dominated by the host galaxy. 

Both values are much larger than the median DM), of 186 pc cm ? derived from popula- 


tion models of FRBs (4). The log-normal distribution of DMy st is also shown in Figure |S4| — 


while very high values of DM;,,,; are disfavored, they are allowed within the model. 


S1.9 Burst energetics 


FRB 20220610A was measured with a peak S/N of 29.8 (which by definition has an uncertainty 
of +1) in CRAFT ICS data formed from 22 ASKAP antennas. The sensitivity in this configu- 
ration scales with the number of antennas, Nant, as VN ant compared to single-antenna searches 
using ASKAP (8), which have a noise level of 2.68 + 0.02 Jy ms for a 1 ms burst (7, 77). We 
have conducted tests on pulsars by comparing ICS to single-antenna mode, and have confirmed 
this scaling to within 97.5 + 2.5%. 

The effective width of a burst in ICS data is increased from its incident (i.e. intrinsic 
plus scatter-broadening) width by DM smearing over the 1 MHz channel bandwidth and by 


the time resolution of 1.182ms. Thus S/N was maximized when the burst was integrated 
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over four time samples. This corresponds to an initial band-averaged fluence estimate of 
S/N x N42? x (At)U5 x Nims = 37Jyms. This is a lower limit on the fluence, some of 
which will lie outside the 4-sample search window used by our real-time detection pipeline. 
Offline fitting of a Gaussian profile to the data produced a full-width half-maximum of 5.6 ms 
and integrated fluence of 43 + 5 Jy ms. Including the potential inefficiency of ICS observations, 
and correcting for the offset from beam center (sensitivity reduced by 2% at central frequency), 


we estimate a band-averaged fluence F, = 45 + 5 Jy ms. 


The isotropic-equivalent energy density of transients is related to the observed fluence via 


An D? 
B (i+ ae (S2) 


The FRB host galaxy has a redshift of 1.016 + 0.002, corresponding to a luminosity distance 
Dy, = 6929 + 2 Mpc using our adopted cosmology (75). We therefore find an energy density of 
E, = (6.4+0.7) x 10? erg Hz~! at the mean emission frequency of 2.543 GHz. 

The spectral behavior of FRBs is complex, with some bursts showing broadband structure, 
and others being isolated in frequency (78). On average, bursts appear to be either brighter (29) 
or more common (4, 79) at lower frequencies. In the former case, fluence scales as F, « v^, 
with a = —1.5*02 (29). If so, then a K-correction of (1+ z)~* = 2.85 should be used to obtain 
E, = (1.8 + 0.2) x 10” erg Hz! at 1271.5 MHz; if not, the K-correction is zero. If the FRB is 
beamed (80), or amplified by gravitational lensing (87), the total energy emitted was less than 


the isotropic equivalent energy. 


S2 Supplementary Text 


S2.1 Nature of the host galaxy 


The spectra in Figure S2]indicate that components (a) and (c) are consistent in redshift to 0.01%, 


—1 


while component (b) is blueshifted by a velocity difference of 150 km s ". With a projected 
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separation of 16 kpc between compact components (a) and (c), they could be two distinct galaxy 
nuclei in a merging or interacting system, while component (b) is the gaseous tail linking them 
and orientated towards the observer. 

The measured (R — K,) colors for the individual components, as well as for the host system 
as a whole (1.70 + 0.18) are bluer than most field galaxies at redshifts around 1 (82—84). At 
z = ] the (R — K,) color approximates rest-frame (U — J) in the AB magnitude system. 
Starburst models (85) only produce (U — J) colors close to these values around 10 Myr after 
the onset of a starburst, or at least 600 Myr later. 

To model the stellar population of the host, we use the stellar population synthesis code 
Prospector (66) which infers the properties of the stellar population and star formation 
history based on spectral energy distribution fitting. We used the integrated photometry from 
the entire complex to make an estimate of the host properties as there is insufficient information 
to characterize individual components. This effectively assumes that all components have the 
same star formation history. 

We fitted the photometry using the stellar population synthesis library python-fsps 
(87, 88). Additional prior assumptions include a Kroupa initial mass function (89), the Kriek 
& Conroy dust-attenuation law (90), mass-metallicity relationship (97), a ratio between dust 
attenuation of old versus young stars, and a continuity model non-parametric star formation 
history (SFH, (92)) using 8 age bins. We used the dynest y (93) package to sample the poste- 
rior probability distribution using a dynamic nested sampling algorithm. The methodology has 
been described elsewhere (94). 

We find the photometry is most consistent with a galaxy with log(stellar mass/solar masses) 
of 9.98* 0-13 at an age of 10.111918 Gyr. Additional parameters are reported in Table 1. We show 
the posterior probability distribution for the stellar mass, marginalized over the other parameters 


in Figure The resulting bimodal distribution could be due to the known degeneracy in 
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spectral energy distribution modeling in age, dust, mass, and metallicity (94). However, given 
the complex morphology of the host it could instead be due to the presence of multiple stellar 


populations corresponding to multiple galaxies. We are unable to distinguish these possibilities. 


$2.2 Implications for the maximum FRB energy 


Determining the maximum energy of FRBs, Emax, requires either a sample of localized bursts, 
or accounting for variations in the Macquart relation. Several parameters — such as source 
evolution, spectral behavior, and the slope of the luminosity function — are expected to be 
correlated with Emax, so must be simultaneously fitted when modelling the data. This approach 
has been taken in previous studies (79, 95, 96), which arrived at mutually consistent results. The 
most precise value reported — derived from a sample of 16 FRBs with measured redshifts, and 
60 without — is log;o (£4. /erg) = 41.26* 027 at 68% confidence (4). This calculation assumed 
an intrinsic emission bandwidth of 1 GHz, and did not include a K-correction, i.e., it allowed the 
rate of FRBs to vary with frequency, but not their energy. That study calculated the isotropic 
equivalent energy, i.e. ignored beaming, which is unconstrained by observations. The burst 
energy function used was an upper incomplete gamma (Schechter) function, with integral slope 
y, and exponential cut-off of the form exp( — £/ E4;), where E is the bandwidth-integrated 
burst energy. In this formalism, the burst energy of FRB 20220610A is (6.4 + 0.7) x 10“ erg, 
3.5 times the cutoff energy — which is allowed but unlikely under that model. 

To test the influence of this observation on the derived Emax, we recalculate the model of (4) 
following a previously described implementation (74, 97), but add FRB 20220610A to their 
sample, and vary Emax only. The resulting likelihoods are shown in Figure|S5| The peak likeli- 
hood is obtained for log; o (Ea /erg) = 41.7+0.2, which is equivalent to log ,)(E,/erg Hz !) = 
32.7+0.2, where the uncertainties were set by assuming a log-uniform prior on Emax to produce 


a Bayesian posterior probability distribution p( Emax). This is 2.7 times higher than the previous 
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estimate, but still within the 90% confidence interval of the previous model. 

Nonetheless, FRB 20220610A is a high-energy FRB. Given the time and frequency resolu- 
tions of the FRB detection algorithm (tsamp = 1.182ms and Av = 1 MHz respectively), the 
DM of 1458.151925 pc cm ^? results in instrumental broadening of the signal over a time dura- 


tion wpm of 5.8 ms at band centre. Including the scattering time 7 and intrinsic width w, the 


total effective width weg becomes 


w? + 7? = 5.9 ms. (S3) 
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Without the wpm, the width would be 1.35 ms. Because observed signal-to-noise scales with 
Wop, the greater width caused by the high DM of the burst decreased our sensitivity by 52%. 

Using the observed properties of FRB 20220610A (DM, width, fluence), the observation pa- 
rameters at time of detection (location in beam, observation frequency, sampling time, number 
of antennas), and estimates of the FRB population parameters, we estimate the probability of 
observing FRB 20220610A at a given redshift, p(z|FRB 20220610A). The results are plotted 
in Figure|S6] We use the best-fitting population parameters (Emax, y, spectral rate scaling index 
a, scaling with the star-formation rate nsf, and log-mean and standard deviation of the host 
galaxy DM jii; and Chost) from the previous model for Ho= 67.4 km s^! Mpc! (75), then set 
Emax to our revised value. We also explore parameter sets corresponding to these same parame- 
ters when at their upper and lower 90% confidence limits. The resulting cumulative probability 
distributions cdf (z|F RB 20220610A) are plotted in Figure|S6| 

Our model for cdf (z|FRB 20220610A) is based on previous work (3, 28), and includes the 
increased probability of encountering halos with increasing z, and the resulting selection effects 
due to DM smearing. However, the model does not include correlations between DM and the 
scattering time 7, because there is no observational evidence for such a correlation (28). It has 


been suggested (98) that a large observed scattering time 7 correlates with a larger value of 
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DMhost- In such a case, a high value of 7 would be evidence for a higher value of DMhpost, 
hence a lower value of DMeosmic and thus z. We regard this as almost certain, however attempts 
to quantify it have been highly uncertain, with at least one FRB with a large excess DM having 
low 7 (99). Such a wide variation between 7 and DM is also found for Galactic pulsars (700). 


Therefore, we regard the approach of using 7 to estimate DM) 5; as having minimal utility. 


$2.3 Limits on depolarization due to multi-path propagation 


The depolarization f due to a spatially variable rotation measure in a foreground scattering 
screen is (74, 101) 
f^1-exp (—2X 50$) ' (S4) 


where f is the observed polarization fraction, A is the observing wavelength and orm parameter- 

izes the RM complexity. This relation assumes the burst is intrinsically 10046 linearly polarized, 

and that depolarization occurs in a medium that is in the foreground to the emitting plasma. 
We use our measured lower limit on the depolarization of FRB 202206104 to set an upper 


limit on egy using the relationship 


nace 
IRM < ee L =os ran (S5) 


assuming f > 0.96. In the frame of the host galaxy this corresponds to ogm < 1.2 rad m ?. 


We regard our assumptions about the depolarization as conservative. If the emission is 
intrinsically less than 100% linearly polarized we would place a smaller upper limit on opm. 
While it is theoretically possible that there is depolarization within the emission region, we 
regard this as unlikely for two reasons. Firstly, there would need to be non-relativistic plasma 
within an emission region comprising relativistic electrons (2, 29). Secondly there would need 
to be rotation measure variations in a region of size S cWo.5msl’ œ% DU x 100km, where c is 


the speed of light, I’ is the unknown relativistic beaming fraction and wo,5 ms is the pulse width 
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in units of 0.5 ms. This would require a very disordered magnetic field, which we regard as 


unlikely. 


S2.4 Potential origins of the excess dispersion 


When considering possible origins for the excess dispersion of FRB 20220610A, we make use 
of additional information encoded by the intervening plasma on the observed temporal scatter- 
ing and rotation measure. In the rest frame of the host galaxy the pulse broadening is 1.88 ms, 
and the rotation measure is 860 + 8rad m^ ?. Because the burst occurred at high Galactic lati- 
tude, the Milky Way contribution to both scattering and rotation measure is small, so we neglect 
it. 

We first consider the possibility that the excess arises in plasma closely associated with 
the burst source (1.e., a circumburst medium). Such an environment has been considered for 
some repeating FRB sources (6, 102). The excess dispersion seen in FRB 20220610A and 
the repeating FRB 20190520B are similar, but the pulse broadening is a factor of 6 higher 
for FRB 20190520B, while the rotation measure of FRB 20190520B is highly variable but 
reaches values more than an order of magnitude above that of FRB 20220610A (103). The 
frequency-dependent polarization seen in FRB 20190520B is consistent with multi-path prop- 
agation through a highly magnetized and turbulent local medium that could explain both the 
scattering and rotation measure; in this model, the rest-frame Faraday dispersion for this source 
is onu = 337radm~?, over two orders of magnitude higher than the limit we derive for 
FRB 20220610A above. Thus, for the local environment to contribute the bulk of the excess 
DM seen in FRB 20220610A, the medium would need to be less turbulent and less magnetized 
than for FRB 20190520B. 

Another possibility is that this excess arises in the interstellar medium (ISM) of the host 


galaxy. We compare the properties of FRB 20220610A in the host galaxy rest frame to those of 
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FRB 20190520B and pulsars in the Milky Way. In Figure[S7] we show the distribution of DM 
and pulse broadening times for pulsars (57, /00) and these two localized FRBs. In the case of the 
FRBs we plot them using the most likely DMyj,, contribution from our DM analysis, given in 
Table [51] FRB 20220610A has substantially lower scattering than expected from a Milky Way- 
like ISM. We have decreased the pulse broadening times of the FRBs by a factor of three to 
account for the plane-wave scattering geometry (104) assumed for FRB pulses scattered in their 
host galaxy, which contrasts to the spherical wave scattering geometry assumed for pulsar pulses 
scattered in the Milky Way. In Figure[58] we show the distribution of DM and rotation measures 
for pulsars (705) and the two FRBs. FRB 20220610A shows a rotation measure similar to most 
Milky Way pulsars at similar DM, but much lower than the RM associated with pulsars close to 
the Galactic Centre (106), some of which have DMs at the upper end of the allowable range for 
FRB 20220610A. 

While the Faraday rotation of the burst is not inconsistent with the relationship between 
dispersion and rotation measure observed in Galactic pulsars, for a Milky Way-like sightline the 
observed excess DM would imply scatter broadening of > 100 ms at 1.4 GHz (Figure[S7). We 
therefore disfavor the host galaxy ISM as the sole source of the excess DM for FRB 20220610A. 

A final possibility is that the excess DM is associated with plasma in an intervening structure 
along the line of sight, but outside of the host galaxy complex. However, we do not expect 
plasma in the intergalactic medium to have high levels of turbulence or magnetization, so is 
unlikely to cause the observed rotation measure or pulse broadening (107, 108). Therefore we 
also disfavor this option as the sole source of excess DM for FRB 20220610A, unless there is 
an intervening galaxy halo, where density and turbulence can be higher. 

Given these constraints, we cannot identify a dominant location for the excess DM seen in 
FRB 20220610A, but the observed properties can be plausibly reproduced by a combination of 


two (or, indeed, all three) of the potential locations. 
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Figure S1: Pulse-averaged spectrum of FRB 20220610A. (A) Spectra for total intensity 7 
(black), Stokes Q (red), Stokes U (magenta), and Stokes V (blue). The variation in flux density 
between Q and U is due to Faraday rotation. The absence of FRB emission above ~1195 MHz 
is due to buffered voltages above this frequency (see text). (B) Linear polarization position 
angle V (black points, with lo error bars). The red line is a maximum-posteriori model, which 


corresponds to a rotation measure of 215 + 2 rad m ?. 
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Figure S2:  VLT/X-shooter spectra of the host galaxy of FRB 20220610A. The Figure is 
the same as Fig 2E, and covers the region of the [O r1] lines with rest wavelengths 3726 and 
3729 A, but for the three components separately and on a velocity scale. The peak of the [O 11] 
3729 À line in component (c) is marked as a vertical dashed line. The red dashed line shows the 
uncertainty in each spectrum. The blue horizontal line shows zero relative flux on an arbitrary 
flux scale. 
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Figure S3: Posterior probability distribution of the host galaxy stellar mass as inferred 
from spectral energy distribution fitting. The posterior probability has been marginalized 
over other parameters. The solid line is the median of the distribution while the dashed lines 
mark the 68% confidence intervals. The distribution shows bimodality that could be due to the 
presence of multiple galaxies. 
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Figure S4: Probability distributions for the rest-frame host galaxy dispersion measures for 
two FRBs. FRB 202206104 is shown as a blue solid line and FRB 20190520B as an orange 
dashed line. The green dotted line shows the best fitting log-normal distribution for DMpost 


derived from a population model (4). 
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Figure S5: Estimation of the maximum FRB energy, Emax for FRB 20220610A. Likelihoods 
L were evaluated after adding FRB 20220610A to the previous sample of (4) and varying Emax 
only. The posterior probability distribution p( Emax) assumes a log-uniform prior on Fmax- 
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Figure 56: Cumulative probability distribution for the redshift, z of FRB 20220610A. This 
was estimated using a previous model (4) (blue solid line); population parameters varied within 
their 90% confidence limits (grey shading); and with our revised estimate of Emax (orange 
dashed line). 
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Figure 57: Dispersion measure and scatter broadening times for Galactic and extragalactic 
sources. Milky Way pulsars are labeled PSR (black) (700, 109), and FRBs FRB 20220610A 
(red) and FRB 20190520B (blue). The error bars on the FRB DM points reflect uncertainties in 
the host DM (see text). Error bars for pulsar DM are too small to display. The solid and dashed 
lines show a mean and 1c relationship derived for pulsars (700). 
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Figure S8: FRB and pulsar dispersion measures and rotation measures. Pulsars (/05) 
close to the Galactic Center (GC) are highlighted in cyan. Other pulsars (labeled PSR) are in 
black (705), while other FRBs (73) are in orange. 
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Table S1: Observed scattering properties of two FRBs inferred to have excess DM: 
FRB 20220610A and FRB 20190520B (6). We scale scattering timescale 7 to 1.4 GHz in 
the host galaxy rest frame, assuming 7 x v ^, to calculate T! ,(444,- The expected contribution 
from the Milky Way, Tmw, is taken from the NE2001 model (57), scaled to v,. We also include 
estimates of DMj,,. from Equation (S1). with 6896 confidence intervals. 


FRB 20220610A FRB 20190520B 


T (ms) 0.511 + 0.012 10+2 
v, (GHz) 1.1475 1.25 
z 1.016+0.002 0.241 0.001 
T! acu (ms) 1.88 + 0.02 12+2 
Tmw (ms) 4.83 x 107 1.87 x 1074 
DMhost (pc cm~?) 13042453 Ip 
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Table S2: Photometric measurements, including the positions and sizes of elliptical aper- 
tures as shown in Figure 2, for each of the components in the host of FRB 20220610A. 


Component: (a) (b) (c) Full aperture 

Right Ascension (J2000) 23^24"17.55* 23h^24"17.645 23524™17.758 23h24"17.65* 

Declination (J2000) —33^30'49.9"  —33?30'48.2" | —33?30'49.9" | —33?30'49.0" 
Position angle (E of N) 45° —35° —80° 0° 
Semi-major axis T 1.6" 0.8" 3” 
Semi-minor axis p 1.1” 0.7” 3” 

g 25.36 +0.06 25.13 +0.06 26.01 0.08 24.15 0.07 

R 24.98 + 0.05 24.82 40.05 2569+0.07 23.78 + 0.06 

J. 23.07 +0.08 23.40 +0.13 23.84+0.12 21.97 X 0.07 

K; 22.92 + 0.09 > 23.66 23.75 +0.14 22.08 + 0.12 

(R- Ks) 2.06 + 0.10 < 1.16 1.94 + 0.16 1.7 + 0.13 
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